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Stable vortex loops in two-species BECs 
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We consider the creation of stable, stationary closed vortex loops, analogue to the vortons and 
superconducting cosmic strings, in cold atom EEC's. We explore the parameter region where these 
solutions are likely to exist and comment on methods to create them experimentally. 



INTRODUCTION 



The existence of vortex solutions is one of the hallmarks 
of superfluidity. They are seen in superfluid ^He as well as 
^ in Bose-Einstein condensates in cold atom traps [1]. Super- 
(30fuid vortices are also the main actors in the leading theory 
^ ' , of rotation and "glitches" observed in rotating neutron stars 
(pulsars) [2] ; superfluid vortices are supported in both the 
Ci hadron and quark matter phases at high density. Similar 
^ structures are also likely to appear in Grand Unified Theo- 
^^ries (GUT) of physics beyond the standard model of particle 
physics. In that context vortex filaments in the early universe 
might serve as seeds around which galaxies could form [3]. 
I The stability of superfluid vortices is guaranteed by their 
topological properties. The velocity circulation around a vor- 
Ch tex is characterized by a phase which wraps axially around 
O the vortex. Because the field must be single- valued this phase 
I ^ I must go through an integer number of cycles and thus cannot 
change continuously. Discontinuous changes through quan- 
tum tunneling or thermal activation are usually negligible. 

Vortices with opposite circulation can annihilate one an- 
other. This phenomenon seems to preclude the possibility 
ly-^ of stable closed vortex loops, because the opposite sides of a 
■^j- closed loop form a vortex/antivortex pair that can annihilate 
* and, in fact, the energetics of the system favors annihilation. 
I Because a vortex has a tension T along its length, the energy 
of a vortex loop is proportional to its length: E ^ 2'kRT, 
T-H where R is the radius of the loop. This energy is reduced 
1[ as the loop shrinks. When the opposite sides of the loop 
. ^ come closer than their thickness they can annihilate. Single- 
species vortex loops can only be stabilized by the Magnus 
J-H force, which will sustain a loop if it is moving at a specific 
velocity in a direction perpendicular to its plane. 

It has been pointed out in the context of relativistic models 
that a stable closed vortex loop can exist if two species are 
competing to condense [1]. Qualitatively, the mechanism 
is straightforward. Suppose in some region there is a bulk 
where the first species was condensed, and there is a vortex 
in that medium. In that vortex's core the condensate of 
species 1 vanishes. If the two species repel each other it will 
be energetically favorable for the second species to condense 
in the core of the vortices of species 1 rather than occupy 




FIG. 1. The arrows denote the momentum (change in phase) of 
species 1 and 2. 



the same space as the first species. Suppose that the vortex 
of species 1 forms a closed loop and, in addition, there is 
a non-zero vorticity of species 2 along the core of the loop 
(see Fig. [T]). The energy cost related to the vorticity I of 
species 2 is proportional to 2ttR{1/R)^ ^ P/R expresses the 
kinetic energy of species 2 inside the vortex. Because this 
contribution increases as the loop shrinks, there will always 
be some nonzero radius where the total energy E ^ TR + 
P / R has a minimum. Thus, if this equilibrium radius is big 
enough to preclude annihilation, the second species provides 
a stabilization mechanism for vortex loops. 

If the particle condensing in the interior of the vortex is 
charged, the vortex will be superconducting and solutions of 
this kind are known as superconducting strings [5] in or- 
der to distinguish them from superconducting flux tubes 
(Abrikosov vortices) where the superconducting region is 
outside of the vortex. Closed loops of superfluid vortices 
stabilized by the mechanism sketched above are know as 
"vortons". Vortons can arise in GUT-scale models during 
the early universe[B]. Quark matter, which can exist in the 
core of neutron stars, also provides the ingredients for the 
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existence of vortons, the role of the two competing species 
played by quasiparticle excitations of the color supercon- 
ducting ground state with quantum numbers of neutral and 
charged kaons [THTU]. 

The purpose of this paper is to study the possibility 
that stable vortons can exist in cold atomic traps. The 
(meta)stability of this kind of solution is analyzed in the 
first part through a calculation of the energy of the vorton 
as a function of its radius and thickness. We first show that 
in some circumstances, the equilibrium radius of the vorton 
is on the scale of microns and is reasonably larger than its 
thickness. We then comment on possible ways of actually 
creating vortons experimentally. 



II. STABILITY AND EQUILIBRIUM PROPERTIES 
OF VORTONS 



Consider a system of two species of spinless bosons. In 
the dilute limit, where the interparticle distances are typi- 
cally larger than the scattering lengths, this system can be 
described by the Hamiltonian 
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Here, (pi and 4)2 are the second quantized fields annihilat- 
ing species 1 and 2, and Mi are the scattering lengths 
and masses of the two species (i=l,2), ai2 is the interspecies 
scattering length, and A/12 is the reduced mass of the two 
species. The chemical potential of species 1 is given, in the 
absence of species 2 and at leading order in the diluteness 
expansion, by 



STrS^riifli 
Ml 



(2.3) 



where, ni is the asymptotic density of the first species. We 
will always be interested in situations where species 1 forms 
the bulk material, so we now eliminate /ii from the dis- 
cussion. We do not include a chemical potential term for 
species 2; instead we will work with a fixed number of par- 
ticles of that type. We will be interested in the phase sep- 
aration regime 4aia2/MiM2 < (ai2/Mi2)^, where the inter- 
species repulsion encourages the two species to stay at sepa- 
rate points in space 

Bose-Einstein condensation is described by a non- 



vanishing matrix element the field operators 0i.2. At low 

— 1/3 

densities [n-^ 2 ^ 01,^2, 012) the mean field approximation 
is valid and the matrix element (also denoted by (/>i,2) satisfy 
the classical equations of motion, that is, the Gross-Pitaevski 
equations. For a straight vortex of only the first species, (pi 
is of the form 



4i{r,e,z) = fi{r)e'' 



(2.4) 



where r, 9^ and z are cylindrical coordinates, /i(f = 0) = 
and fi(r ^ (X)) ~ \fni. To keep (pi single valued the winding 
number j must be an integer. We always use j = 1, as 
vortices with larger j are unstable and split into j vortices 
with winding number 1, but keeping general j is useful for 
studying the energy of a vorton. The precise profile /i(r) 
can be obtained by solving the Gross-Pitaevski equations or, 
equivalently, minimizing the energy. 

Simple scaling arguments imply that the solution /i(r) 
will change from /i = to /i = \/ni over a distance of 
order b ~ l/-\/87rniai, the string thickness. The toroidal 
geometry of the vorton makes an analogous calculation a lit- 
tle involved. We can, however, bypass most of the difficulty 
by assuming that the radius i? of the vorton is much larger 
than its thickness 8. In this case we can compute the energy 
of the vorton by computing the energy per length (that is, 
the tension) of a straight vortex and multiplying it by the 
length 27ri?. This approximation neglects the energy associ- 
ated with the curvature of the vorton. For a straight vortex 
of the first species with an internal current of the second, the 
mean-field solution has the form 



^2 = /2(r)e^'=^ 



(2.5) 



To encode the fact that we wish to consider a vortex loop 
of radius i?, we impose periodic boundary conditions in z, 
identifying z = ±nR, so that z spans the arclength of the 
loop. The continuity of the phase of (/)2 implies that k = l/R, 
with integer I. 

Plugging the functional forms in (2.5 1 with k — l/R into 



the Hamiltonian given in (2.1 1, the total energy of a vorton 
of a radius R is 




Mr) 




+ V{flj2) + 



which becomes more accurate as the ratio R/S grows. One 
might worry that including curvature corrections might en- 
courage the vorton's radius to be smaller. Brief considera- 
tion alleviates this fear: the configuration we consider when 
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formulating (2.6) has no gradients of 0i along the vorton's 



length. Curvature effects will make those terms nonzero, and 
the nonuniformity of (f)i will provide a considerable potential 
barrier. Therefore we expect that calculation that included 
curvature effects would find an R larger than the R wc find 
from considering the energy expression in (2.6). 



In (2.6) we measure the energy in relation to the homo- 
geneous ground state j — 0, fi — /2 = 0. These 
parameters correspond to no vorton at all. Any vorton will 
have a greater energy, and thus must be at best metastable. 

We now argue that vortons are indeed long-lived. If a 
vorton is initially created with a non-zero number N2 parti- 
cles of species 2 in its interior, we hope they remain trapped 
there, or the stabilizing energy that scales inversely with R 
will disappear, causing the vorton's collapse. Escaping the 
vorton would require these particles to go through a large 
region where species 1 is condensed, and that is energeti- 
cally expensive. One could also imagine a bubble of species 
2 detaching from the vorton and moving towards the edge of 
the bulk, but that costs an amount of energy proportional 
to the area of the bubble due to its surface tension. These 
escape mechanisms are costly, and thus we treat N2 as a con- 
served quantity that is confined to the vorton's interior. To 
implement this conservation, when we minimize the energy, 
we keep 



N2 



2TrR j dr 27rr|/2(r)|^ 



(2.7) 



fixed. 

To summarize, we find our vorton candidates by numeri- 
cally minimizing the energy in (2.6 1 in relation to fi{r), f2{f) 
and R while keeping the parameters ni, ai, 02, ai2. Mi 
and M2 fixed. We have eliminated /ii, and henceforth take 
j — 1. A typical example of the profile functions /i(r) and 
/2 ('') resulting from this minimization is shown in Fig. [2j 

It is easy to foresee some qualitative dependences before 
doing the numerical minimization. For instance, the larger 
the species 2 vorticity I is, the larger R and R/S should 
be. A larger value of 012 increases interspecies repulsion and 
shrinks the region where both species coexist. This tendency 
should quickly saturate, as the coexistence region becomes 
negligible and 012 becomes effectively infinite. We can also 
see that if we increase ni by a factor of rj and reduce the 
number of particles of species 2 by a factor of ,yf], both the 
radius and thickness decrease by a factor of ^yf| and therefore 
R/S remains unchanged. 

On the other hand, the dependence of R/S on some pa- 
rameters are harder to anticipate. For example, as M1/M2 
grows, the terms in the energy describing species 2 (and its 
interaction with species 1) grow. Because species 2 exerts 
a pressure counterbalanced by the bulk, this increase of the 
contribution from species 2 should increase both R and S, 
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FIG. 2. The radial density profile of the two species ^Li and ®^Rb 
and I = 5 The remaining parameters are given in the first row of 
Table E 



so it is difficult to know in advance if their ratio will grow 
or shrink. Numerical studies indicate that the ratio R/S in- 
creases with Mi/M2- We also observe an increase of R/5 
with increasing oi and decreasing 02. 

To make our detailed discussion more concrete, we study 
the vortons formed when the roles of species I and 2 are 
played by the \F = l^mp = 1) hyperfine state of ^^Rb 
and |1,I) state of ^Li, respectively. The interspecies scat- 
tering length can be tuned through the use of Feshbach res- 
onance p[^. A naive estimate of thickness of the vorton 5 
is simply thickness of a *''Rb vortex, which is given by the 
healing length 1/ ^/Etttv^oi, an order-of-magnitude estimate 
supported by our numerical calculations. For definiteness, 
we pick the parameters which correspond to the first row of 
Table |l] and estimate the healing length to be 0.45 microns, 
while S = SAfim. The ratio of these two quantities is roughly 
7.5, which is seen in Fig. [2] 

We can estimate N2 by simply multiplying the cross- 
section of the vorton, its length, and the density of the 
lithium, N2 ~ (7r(5^)(27ri?)n2. As an example, for Lithium 
and Rubidium, if the density of lithium is roughly the same 
as the density of rubidium (10^'^/cc), and we want R of the 
order of 10 microns with an R/S of 2.5, then N2 should be 
on the order of 10'*. We show some example situations in 
Table ID where we list the results for Li-Rb {M1/M2 = 12) 
vorton and K-Rb {M1/M2 = 2.12) vorton for some realis- 
tic values of ni and N2. These results indicate that the 
vorton has only a moderate sensitivity to the number of par- 
ticles present. It also indicates that the R/S ratio is not 
large, making the corrections to the approximation we used 
in computing its energy sizable. However, as argued, the cur- 
vature effects increase R/S, so the existence of stable vortons 
is nonetheless supported by this analysis. 
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TABLE I. R/S for diflerent values of parameters for I = 5. The 
values of scattering lengths used were obtained from [13) . |14| . 
|15| . The scattering lengths are measured relative to the Bohr 
radius ao. 



III. PRODUCTION MECHANISMS 

Up to now we have discussed the properties of static vor- 
tons. We now briefly comment on two possible ways to ac- 
tually create them experimentally. 

Raman Scattering and Gauss-Laguerre Beams A possible 
way to engineer a vorton in cold atoms traps suggested by 
J.V. Porto[16j involves three ingredients. 

First, a closed vortex loop of *^Rb can be created by the 
use of two counterpropagating light beams, each one with a 
central circular region with one frequency and an outer annu- 
lar region with another frequency, as shown in Fig. [3j In the 
left-moving beam the inner frequency uja is chosen to match 
the Di line of rubidium. The outer frequency wt is chosen 
to match a transition from the same excited state to another 
hyperfine state of *^Rb. In the counterpropagating beam, 
the frequencies Wa and LOb are exchanged. Thus rubidium 
atoms located in the inner region of the beams will absorb 
a photon with momentum hk = huJa/c and emit a photon in 
the opposite direction with nearly the same momentum, due 
to the stimulation of the second beam. The net effect is that 
the ^^Rb atoms in the inner region acquire a momentum 2hk. 
In the annular outer region of the beams the emission and 
absorption are reversed and the atoms acquire a momentum 
-2hk. 

At the boundary separating the two regions the shearing 
of the *^Rb encourages vortex loop formation. The recoil 
energy of the rubidium atom, for optical photons, is of the 
order of 10 ^lK so heating may be an issue. This problem 
might be circumvented by having only a small fraction of 
the atoms go through the absorption-emission process. Fur- 
ther scattering of these atoms will distribute its energy and 
momentum among nearby atoms. 

The second ingredient is that the creation of the rubidium 
loop be accompanied by a change in magnetic field to a value 
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FIG. 3. Geometry of the two countermoving beams. The arrows 
denote the direction the beams are moving and the letters "a" 
and "b" stand for the two optical frequencies. 



close to the Feshbach resonance between the rubidium and 
lithium atoms. Such a tuning will entice the system to phase 
separate, so that the lithium will seek locations with little 
rubidium: either the boundary of the bulk or the interior of 
the rubidium vortices. 

Finally, the lithium atoms must acquire a net angular mo- 
mentum around the vorton before it has time to collapse. 
This can be accomplished with the use of a Gauss-Laguerre 
beam. The small values of / < 5 considered in the numeri- 
cal examples in Table |l] were motivated by the limitations of 
the Gauss-Laguerre beam technique. With improved Gauss- 
Laguerre techniques, a higher I and larger R/5 could be 
achieved. 

The Kibble-Zurek Mechanism At high temperatures the 
phase of the bosonic fields 0i and 02 at different points in 
space is uncorrelated. At low temperatures, in the super- 
fluid state, the phase has long range correlations. However, 
following a rapid quench, the phases do not have enough 
time to correlate and defects form at the regions where dif- 
ferent phases meet. This process (the Kibble-Zurek mech- 
anism [T71 US] ) was demonstrated to create regular vortices 
in ^^Rb[T9^. In the present case two ingredients should be 
added. Vorticity needs to be given to the ""Li using standard 
methods and a change of the external magnetic field to the 
Feshbach resonance, in order to impose phase separation. 

Any rubidium vortex generated by the standard Kibble- 
Zurek mechanism will be filled with lithium atoms and any 
vorticity of the lithium enclosed by the vortex loop will be 
conserved and will guarantee the stability of the vorton. 

Some simplifying assumptions allows us to estimate the 
probability of forming a vorton during such a quench. First, 
we assume that the quench is sufficiently rapid such that 
the correlation length of the phases of each field is the same 
as before the quench. Second, we assume that there is lit- 
tle interaction between the two species before the quench 
(which can be enforced by tuning the interspecies scattering 
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length to zero) but, after the quench, lithium atoms imme- 
diately fall into the potential well generated by the core of 
rubidium vortices. The typical distance between rubidium 
vortices after the quench is given by the correlation length 
^1 ^ h/k ~ MiT / ^Ai^aini, where k is the typical momen- 
tum of a phonon, related to the typical energy of a phonon 
e ^ T ^ ck, where c = V^^nxini/Mi is the speed of sound. 

Since the orientation of the rubidium vortices is random, a 
fraction of order one of those vortices will connect to a vortex 
with the opposite orientation and form a vorton. In those 
cases, lithium atoms will be trapped in the core of the closed 
vortex and guarantee its stability as long as there is vorticity 
in the helium along the closed vortex loop. We can estimate 
this vorticity by a standard argument |18j. The vorticity I is 
the integral of the phase of 02 along the vorton line 

I = j ds.V92, (3.1) 

where 02 — ^26*^^- The phase of 02 changes over a distance 
of order ^2 ^ fi/k ^ M2T / ^/4mi2n2, in exact analogy with 
01. If we assume that the phase of 02 is uncorrelated in 
different regions apart by more than ^2 , we expect the phase 
of 02 to change between and tt a number of times simply 
given by 2ttR/£^2- If the phase always changed in the same 
way, that would be our estimate for I. However, since the 
phase can change either clockwise or counterclockwise, we 
expect / to scale like a one-dimensional random walk, 



1 = J ds.V02-y^, (3.2) 

where R is the vorton radius. The radius of the vorton is of 
order R so we estimate 

^-\i^-\f^( — r- (3-3) 
V 6 V M2 \ainj 

This estimate suggests that we should pick atoms with very 
different masses, and that for increased stability created by 
the Kibble-Zurek mechanism, one should begin with the light 
species at least as dense as the heavy species, but more dense 



if possible. We find that this heuristic estimate suggests, for 
the case of an equal density Rb/Li mixture and scattering 
lengths given by the first line of Table |lj vorticity of Z « 7, a 
value similar to I = 5 which was considered in our stability 
analysis. 

IV. SUMMARY AND OUTLOOK 

We considered a cold atom instantiation of vortons, two- 
species vortex loops. While vortex loops of a single species 
will collapse unless supported by the Magnus force, is is pos- 
sible for vortons to achieve a stable equilibrium radius at rest 
with respect to the bulk. Vortons are supported by several 
astrophysical systems, but have never been observed experi- 
mentally in any context. For realistic laboratory conditions 
it should be possible to create vortons with radii on the order 
of micrometers, with R/S large enough to provide a believ- 
able toroidal geometry. 

Aside from the relations to astrophysical systems, such 
solitons are an interesting aspect of two-species BECs in their 
own right. It is an appealing notion to think that such rich 
structure might be created and manipulated in the labora- 
tory. The requisite two-species BEC cold-atom traps are a 
well-established technology, and many two-species pairs have 
been trapped (see, eg. pOH22]). While we have focused on 
a lithium-rubidium system for numerical results, in princi- 
ple any two-species BEC supports such structures for some 
choice of vorticity and initial densities. Aside from comput- 
ing equilibrium properties, we have discussed two possible 
mechanisms for crafting vortons in cold atom traps. Creat- 
ing such objects would demonstrate their feasibility in other 
contexts, provide an exciting playground for further investi- 
gation, and would widen the possible collaborative avenues 
between astrophysics and cold atom physics. 
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